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The effects of flame structure on the extinction limits of CH 4 -0 2 -N 2 counterflow diffusion flames were 
investigated experimentally and numerically by varying the stoichiometric mixture fraction Z st . Z u was 
varied by varying free-stream concentrations, while the adiabatic flame temperature was held fixed by 
maintaining a fixed amount of nitrogen at the flame. Z st was varied between 0.055 (methane-air flame) and 
0.78 (cbluted- methane-oxygen flame). The experimental results yielded an extinction strain rate K ^ of 375 
s - 1 for the methane-air flame, increasing monotonically to 1042 s " 1 for the diluted-methane-oxygen flame. 
Numerical results with a 58-step Cl mechanism yielded 494 s' 1 and 1488 s’ 1 , respectively. The increase 
in K nt with Z sl for a fixed T ad is explained by the shift in the 0 2 profile toward the region of maximum 
temperature and the subsequent increase in rates for chain-branching reactions. The flame temperature 
at extinction reached a minimum at Z rt » 0.65, where it was 200°C lower than that of the methane-air 
flame. This significant increase in resistance to extinction is seen to correspond to the condition in which 
the OH and O production zones are centered on the location of maximum temperature. 


Introduction 

Studies of strained laminar diffusion flames have 
proven useful, both for developing a fundamental 
understanding of flame structure and extinction and 
for application in the flamelet model of turbulent 
combustion. Flame extinction resulting from exces- 
sive strain is of particular importance and has re- 
ceived considerable attention. For a strained diffu- 
sion flame, Lirians [1] results for one-step 
irreversible kinetics have shown that, for a given 
free-stream temperature, the oxidizer or fuel can 
leak through the flame. Which reactant leaks 
through depends on the ratio of oxidizer to fuel mass 
fractions. In the present work, the stoichiometric 
mixture fraction, 

Z st = (1 + Y p^'WoVo/y'o, - f) 1 (1) 

is used to characterize this ratio, where Y is the mass 
fraction, W is the molecular weight, v is the stoichio- 
metric coefficient, and the subscripts O and F refer 
to oxidizer and fuel, respectively. For Z st < 0.5, ox- 
idizer leaks through, whereas for Z st >0.5, fuel leaks 
through. A different reactant penetrates the flame 
with (Afferent Z st because with one-step kinetics, the 
controlling feature for reactant leakage is the global 
structure of the flame, which is dictated only by the 
composition of the free-streams (assuming die tem- 
peratures of the free-streams are equal). 

For the methane-air flame (Z st = 0.055), it is well- 
known that oxygen leaks through the flame while 
methane is completely consumed [2]. One-step 


asymptotic analysis incorrectly predicts fuel leakage 
because, although the basic flame structure suggests 
fuel leakage, the rate-limiting oxidation kinetics re- 
sult in oxygen leakage. Reduced mechanisms with 
two steps or more can capture separately the kinetics 
of fuel and oxygen consumption and thus predict 
oxygen leakage [2]. Although Lifian s results do not 
correctly identify which reactant leaks through the 
flame, they do demonstrate that extinction is inti- 
mately linked with outer flame structure, as dictated 
by Z st , such that when the basic structure of the 
flame is altered, the extinction strain rate will be af- 
fected. 

A common example of the effects of free-stream 
composition on extinction is the decrease in that 
occurs when reactants are diluted [3-5]. Fuel dilu- 
tion decreases the flame temperature but also in- 
creases the stoichiometric mixture fraction, moving 
the flame away from the oxidizer and toward the 
fuel. The reduction in flame temperature slows 
chain branching while only weakly affecting radical 
recombination, making it easier to extinguish the 
flame. However, the basic structure also changes due 
to the shift in the flame location. To the best of the 
authors’ knowledge, a systematic study of the effects 
of the shift in flame location on flame structure and 
extinction has not been performed. 

Ishizuka and Tsuji [4] studied the effect of inert 
addition on extinction of methane flames. They di- 
luted either the fuel or the air stream to determine 
the limiting fuel and oxidizer concentrations for sta- 
ble burning. Their results showed that there is a 
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TABLE 1 

Extinction Conditions for CH 4 -0 2 -N 2 Flames with T ad = 2226 K. 


Ych< at Z = 1 
Y 0a at Z = 0 
Kert (s- 1 ) (exp.) 
(s' 1 ) (num.) 

T ert (K) 


Flame A 


0.055 
1.0 
0.233 
375 ± 10 
494 
1770 


Flame B 


0.35 
0.157 
0.339 
482 ± 10 
814 
1752 


Flame C 


0.50 
0.111 
0.437 
601 ± 10 
996 
1651 


Flame D 


0.65 
0.0845 
0.632 
770 ± 20 
1226 
1572 


Flame E 


0.78 

0.0709 

1.0 

1042 ± 20 
1488 
1656 


limiting flame temperature that constrains the min- 
imum concentrations for the reactants. Puri and Ses- 
hadri [3] studied the effects of dilution on propane 
as well as on methane flames. In this study, Z st was 
held constant as oxygen mass fraction was varied to 
identify the extinction strain rate. A small range of 
Z st was considered. Chen and Sohrab [5] also studied 
the limiting concentrations for methane flames and 
extended the study to include pure oxygen in the 
oxidizer, thus obtaining-limiting conditions under a 
wide range of fuel and oxidizer concentrations. In 
all of the mentioned studies, the free-stream com- 
positions were varied with no attempt to maintain 
constant flame temperature. Thus, the extinction 
conditions observed are a consequence of changes 
in flame temperature and global structure. 

The objective of this study is to identify the spe- 
cific influence of flame structure (flame location) on 
extinction, and to this end, Z sl is varied while main- 
taining a constant adiabatic flame temperature. The 
extinction strain rates are obtained both experimen- 
tally and numerically. The experimental extinction 
strain rates serve to validate the numerical results, 
and the numerical results are then used to identify 
the structure and to explain the effects of Z st on ex- 
tinction. 

This approach for studying the effects of Z st at 
constant flame temperature has been employed by 
Du and Axelbaum [6] to understand the effects of 
structure on soot inception in ethylene flames. Lin 
and Faeth [7] have since employed this approach to 
experimentally study soot inception and flame ex- 
tinction for a variety of fuels, and Sung et al. [8] 
followed this approach to understand the endother- 
mic reactions on the fuel side of methane diffusion 
flames. 


Experimental and Numerical Methodology 

The experimental apparatus is described in detail 
in Du and Axelbaum [6]. Briefly, a flame is estab- 
lished between two 11-mm opposed jets spaced 8 
mm apart. The jets emanate from tubes with 
honeycomb cores placed 50 cm upstream of each 


exit. An annular co-flow of nitrogen is added to both 
fuel and oxidizer to eliminate oxygen entrainment 
from ambient air and to extinguish the flame that 
would exist outside the region of interest. The com- 
plete flow field has not been characterized, but it is 
expected to be between potential and plug flow [9]. 

To obtain extinction strain rates, the flow rate is 
increased until extinction occurs. A bypass is used to 
ensure that the composition remains fixed as strain 
rate is increased. The flame is then reestablished at 
conditions very near extinction, and the velocity pro- 
file along the stagnation streamline is measured with 
laser-Doppler velocimetry. The strain rate, K t is 
identified as the magnitude of the slope of the ve- 
locity profile upstream of the flame. Because the 
densities of the fuel and oxidizer free-streams are 
similar, the velocity gradients on both sides of the 
flame are identical within experimental uncertain- 
ties. 

The numerical scheme employed was developed 
by Smooke et al. [9]. The version of the code used 
here assumes that the flow field outside the mixing 
layer can be described by potential flow. The reac- 
tion mechanism used is the a 58-step Cl mechanism 
given in Ref. 10. 


Results and Discussion 

To study the effects of flame structure, Z sf must 
be varied while maintaining a constant adiabatic 
flame temperature. Z s( can be varied by simply add- 
ing inert to one of the reactants. However, as noted 
above, this will also change the flame temperature, 
and in this work, the objective is to isolate the effects 
of flame structure resulting from shifts in flame lo- 
cation. Therefore, our approach is to consider flames 
with different values of Z st but with the same adia- 
batic flame temperature T a j. To accomplish this, we 
start with the methane air flame as the reference 
condition. The stoichiometry is 

CH 4 + 20 2 + 7.52 N 2 -4 products (2) 
The adiabatic flame temperature for this reaction is 
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FiC. 1. Numerical results for maximum temperature as 
a function of 1/K for flames A (Z J( = 0.055) and E (Z st — 
0.78). For clarity, only select data are shown near K exi . 



FiC. 2. Extinction strain rate and maximum temperature 
at extinction as functions of stoichiometric mixture frac- 
tion. Open squares denote experimental data for K ext , and 
curves denote computational results. 


2227 K. Z st can be varied without changing T ac i sim- 
ply by maintaining a fixed amount of nitrogen while 
varying the relative amount that is added to the fuel 
and oxidizer. Experimentally, this is accomplished by 
choosing flow rates that satisfy 

Q^ofQo + Qu2.^Qf = 3-76 (3) 

where is the flow rate of nitrogen mixed with 
the ith species, and subscripts O and F refer to 0 2 
and CH 4 , respectively. To obtain the widest range of 
Z st and thus more clearly delineate structural effects, 
flames are considered where all the nitrogen is in- 
troduced with the oxidizer (methane air flame, Z st 
« 0.055) to the other extreme where it is introduced 
with the fuel, such that the oxidizer is pure oxygen 
(Z st - 0.78). 

The validity of this approach to maintaining con- 
stant flame temperature is demonstrated in Fig. 1, 
where the maximum temperature obtained from the 
numerical calculations is plotted as a function of 1/ 
K for Z st = 0.055 and for Z st = 0.78. As 1/K 


that is, in the limit of infinite Damkohler number, 
the temperatures of both flames asymptote to a com- 
mon value, which is as expected because the adia- 
batic flame temperatures are the same. Nonetheless, 
the behaviors near extinction vary considerably. For 
low Z st , the flame temperature is nearly constant un- 
til very close to extinction, where it abruptly falls off. 
For large Z st , the temperature falls off much more 
gradually and persists to a much lower temperature. 
Though the fuel and T a j are the same, the response 
of the flames to increasing K are very different. 

The experimental extinction strain rates, averaged 
over five separate measurements, are listed in Table 
1 for five different flames, flames A-E. For the 
methane air flame, = 375 s -1 , which is in good 
agreement with the results of Chelliah et al. [11], 
wLo found K^t — 380 s" ! . The extinction strain 
rates are seen to increase with Z st , even though T a j 
is constant. For example, the strain rate required to 
extinguish the diluted methane-oxygen flame is al- 
most three times as high as for the methane-air 
flame. 

Table 1 also includes numerical results for and 
the flame temperature at extinction, T^. In Fig. 2, 
experimental and numerical results for are plot- 
ted as functions of Z st . The numerical results are 
found to correctly predict the increase in K^ t with 
Z st , yet they consistently overpredict by 30-50%. 
This discrepancy for the methane air flame has been 
observed by others and is due to the simplified ki- 
netic scheme employed in the model and the as- 
sumption of potential flow outside the mixing layer, 
when the actual flow field is somewhere between 
potential and plug flow [11]. For the purposes of this 
study, a precise modeling of the flow field is not nec- 
essary and the potential flow assumption should be 
adequate. Thus, the qualitative agreement between 
experiment and theory over this wide range of Z sl 
demonstrates that C2 kinetics are not needed to pre- 
dict the observed trends, and that the Cl scheme 
employed is sufficiently valid to allow it to be used 
to extract an understanding of how Z st affects the 
structure of methane flames. 

Before considering the results from detailed ki- 
netics, it is instructive to evaluate how the outer 
structure of the flame is affected by increases in Z st 
by considering a simple flame sheet argument. Since 
the adiabatic flame temperatures of flames A-E are 
the same, the primary effect of increasing Z st is to 
shift the flame from the oxidizer toward the fuel. In 
doing so, the outer structure of the flame is affected. 
Considering the equilibrium (flame sheet) solution 
in Z space, shown in Fig. 3, the changes in outer 
structure are readily interpreted. For purposes of 
discussion, let T ^ represent a temperature that 
would define the boundaries of the reaction zone if 
the kinetics were finite. Noting the temperature pro- 
files and the nominal boundaries of the effective re- 
action zone, it follows that for flame A, there is, on 
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Mixture Fraction, Z 


FlC. 3. Flame sheet solution for temperature and reac- 
tant profiles in Z space for flames A and E. 




Mixture Fraction, Z 


FlC. 4. Numerical results for methane, oxygen, and tem- 
perature profiles in Z space as functions of strain rate K for 
(a) flame A, and (b) flame E. 
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FlC. 5. Numerical results near extinction for major spe- 
cies profiles and temperature in the physical coordinate for 
(a) flame A, and (b) flame E. 

the average, more fuel than oxidizer in the high-tem- 
perature reaction zone for the low Z st flame [2], 

Now, considering flame E, which has the highest 
value of Z st , it is clear that the average oxygen mass 
fraction within the reaction zone is much higher than 
that of flame A. The higher mass fraction within the 
high-temperature region of the flame would be ex- 
pected to accelerate branching chemistry and pro- 
duce more radicals. Thus, the flame would be 
stronger and more difficult to extinguish. 

The validity of the preceding interpretation is con- 
firmed by the numerical investigations employing 
detailed kinetics and transport properties. Figure 4 
shows the reactant and temperature profiles ob- 
tained for flames A and E, plotted in mixture fraction 
space to allow for convenient comparison with Fig. 
3. The high strain rate curves in Fig. 4 are at their 
corresponding extinction limits. These results show 
that oxygen leakage is common to all flames, as is 
the absence of methane leakage. Therefore, the ox- 
ygen consumption reactions are rate limiting at high 
Z st , as they are for the methane air flame. Other than 
the increase in oxygen leakage, the methane and ox- 
ygen profiles are little affected by increasing strain 
rate. The temperature falls off, and the peak tem- 
perature shifts towards smaller Z. 

Figure 5 shows the profiles of temperature and 
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FlG. 6. Numerical results near extinction for profiles of 
OH, O, H, and temperature in the physical coordinate for 
(a) flame A, and (b) flame E. 


major species for conditions near extinction in the 
physical coordinate for flames A and E, with the lo- 
cation of the stagnation point shown for reference. 
As Z st increases, the flame location moves toward the 
stagnation point, with flame C (Z sl = 0.5) being lo- 
cated approximately at the stagnation point. With 
further increases in Z st , the flame location moves 
past the stagnation point toward the fuel boundary. 
In addition, the amount of oxygen in the high-tem- 
perature region increases as Z st increases. For ex- 
ample, at T = 1600 K on the oxidizer side, there is 
5.7 mole % 0 2 for flame A as compared to 23.0 mole 
% for flame E. Thus, these findings are consistent 
with the flame sheet arguments described previ- 
ously. 

The amount of fuel in the high-temperature zone 
is considerably greater in flame A than in flame E. 
Although the rate of fuel pyrolysis is not expected to 
have a dominant effect on flame extinction, this find- 
ing does have important implications for soot for- 
mation; and in Ref. 6, the effects of flame structure 
on soot formation have been discussed. 

In Fig. 6, the profiles of O, H, and OH near ex- 
tinction are shown for flames A and E. The concen- 
trations of OH and O are higher in flames with 
higher Z st> whereas those of H are somewhat lower. 
The combination of high temperature and large 




FiC. 7. Reaction rate profiles for the primary radical pro- 
duction reactions and temperature in the physical coordi- 
nate for (a) flame A, and (b) flame D. 


oxygen concentration in the vicinity of the hydrogen 
radical pool accelerates the primary chain branching 
reaction, 

H + 0 2 -> OH + O (4) 

because this reaction has a high activation energy 
and is first order in 0 2 concentration. Thus, higher 
Z st favors production of OH and O and consumption 
of H. Numerical results for flames with K = 1/2K**, 
have also been obtained to ensure that the given 
trends are not anomalies associated with the unsta- 
ble conditions at extinction. Though the radical pool 
increases when K is decreased, the trends are con- 
sistent with those of Fig. 6 and lead to the same 
conclusions as to the effects of Z st . 

The Z st = 0.78 flame can withstand a lower tem- 
perature before extinction because the increases in 
OH and O provide greater resistance to extinction, 
allowing the flame temperature at extinction to 
be lower. As shown in Fig. 2, the dependence of T ^ 
on Z st is not monotonic. T ^ approaches a minimum 
of 1299°C near Z 5t = 0.65 (flame D), which is 200°C 
lower than for the methane-air flame. Figure 7, 
which is a plot of rates for the radical production 
reactions for flames A and D, suggests an explana- 
tion for this behavior. The reaction rates for these 
reactions, particularly Eq. (4), are centered at the 
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location of maximum temperature for flame D, 
whereas they are far from this location for flame A. 
Thus, the temperature where the hydroxide and ox- 
ygen radicals are produced is comparable to the 
flame temperature for flame D, as opposed to about 
100°C less for flame A. Consequendy, the flame tem- 

S erature can be reduced more for flame D than for 
ame A before the temperature in the radical pro- 
duction zone is reduced to an extent that the flame 
can no longer support combustion. For Z st > 0.65, 
the radical production zone is shifted farther toward 
the fuel, and the maximum temperature must again 
be higher to sustain radical production. 


Conclusions 

For the CH 4 -0 2 -N 2 system at a given adiabatic 
flame temperature, the extinction limits increase 
with stoichiometric mixture fraction. This increase, 
which results from the shift in the 0 2 profile into 
regions of higher temperature, raises the OH and O 
production rates and concentrations in this region, 
yielding stronger flames. This shift in the relative 
temperature and 0 2 profiles also results in a differ- 
ence of up to 200°C in the extinction temperature 
for the same fuel and the same adiabatic flame tem- 
perature. This result is contrary to what is suggested 
from the data of Ishizuka and Tsuji [4], that there is 
a unique limit temperature that constrains the extent 
of dilution that a flame can sustain. The position of 
the radical production zone relative to the maximum 
temperature appears to be critical in identifying the 
limit temperature for stable burning. 

The results of this study are also relevant to prac- 
tical flames wherein the fuel is diluted and/or the 
oxidizer is enriched. There are a number of appli- 
cations where enriched oxygen is presently used, 


notably in blast furnaces. Enriched oxygen combus- 
tion has also been proposed to reduce particulates, 
for example, in diesel engines [12]. 
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